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A naive computation of the correlation functions of fluctuations generated during inflation suffers from log- 
arithmic divergences in the infrared (IR) limit. In this paper, we propose one way to solve this IR divergence 
problem in the single-field inflation model. The key observation is that the variables that are commonly used 
in describing fluctuations are influenced by what we cannot observe. Introducing a new perturbation variable 
which mimics what we actually observe, we propose a new prescription to solve the time evolution of pertur- 
bation in which this leakage of information from the unobservable region of the universe is shut off. We give a 
proof that IR divergences are absent as long as we follow this new scheme. We also show that the secular growth 
of the amplitude of perturbation is also suppressed, at least, unless very higher order perturbation is discussed. 



I. INTRODUCTION 

Inflation has become the leading paradigm to explain the 
seed of inhomogeneities of the universe as seen in the Cos- 
mic Microwave Background (CMB). Despite its attractive as- 
pects, there are still many unknown aspects about inflation 
scenario 01 S HI 01- When we discuss the primordial fluc- 
tuations within linear analysis, many inflation models predict 
almost the same results, which are compatible with the obser- 
vational data, although the underlying models are quiet differ- 
ent. To discriminate between different inflationary models, it 
is important to take into account nonlinear effects |5j 



271 12811 . However, it is widely recognized that we encounter 
divergences originating from the infrared (IR) corrections in 
computing the nonlinear pert urbations generated during infla- 
tion jm n m m sn n m, m m m m. These 

divergences are due to the massless (or quasi massless) fields 
including the inflaton which gives the almost scale invariant 
power spectrum, i.e., T{k) oc k~ 3 . 

We can easily observe the appearance of logarithmic diver- 
gences in the IR limit from the direct computation of loop cor- 
rections under the assumption of scale invariant power spec- 
trum. As a simple example, let us consider a one-loop dia- 
gram containing only one four-point vertex as shown in Fig.Q] 




FIG. 1: One-loop diagram having one four-point interaction vertex 
for the two-point correlation function. 



The end points of the loop are connected to the same four- 
point vertex. Therefore the factor coming from the integral of 
this loop becomes J d 3 kJ'(k). Substituting the scale invariant 
power spectrum into CP(fc) oc k~ 3 , we find that the integral is 
logarithmically divergent in the IR limit like J d 3 k/k 3 . As is 
seen also in this simple ex ample, the IR divergences are typ- 
ically logarithmic J23 IH [ll [M H3, ill IH 53] ■ To 
be a little more precise, we also need to care about UV diver- 
gences. However, since the fluctuation modes whose wave- 
length is well below the horizon scale (sub-horizon modes) 
do not feel the cosmic expansion, they are expected to behave 
as if in Minkowski spacetime. Namely, the quantum state 
of sub-horizon modes is approximately given by the one in 
the adiabatic vacuum. Hence, the sub-horizon modes will not 
give any time-dependent cumulative contribution to the loop 
integral after appropriate renormalization. They are therefore 
irrelevant for the discussions in this paper. Throughout this 
paper, we neglect the contribution due to sub-horizon modes 
by introducing the UV cut-off of momentum at around the co- 
moving horizon scale aH where a is the scale factor and H is 
the expansion rate of the universe. 

As a practical way to make the loop corrections finite, we 
often introduce the IR cut-off at the co-moving scale cor- 
responding to the Hubble horizon scale at the initial time, 
ctiHi 114 111 . This kind of artificial IR cut-off is not fully sat- 
isfactory because it leads to the logarithmic amplification of 
the loop corrections as we push the initial time to the past like 



(Loop integral) 



d 3 kk~ 3 oc \og(a/ai) , (1.1) 

aH>\k\> ai H 
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where <Zj is the scale factor at the initial time and we neglected 
the time dependence of the Hubble parameter. Due to the non- 
vanishing IR contribution, the choice of the IR cut-off affects 
the amplitude of loop corrections. Furthermore, the reason 
why we select a specific IR cut-off is not clear. This means 
that, in order to obtain a reliable estimate for the IR correc- 
tions, we need to derive a scheme to make the corrections fi- 
nite from physically reasonable requirements. This is what we 
wish to discuss in this paper. 

To begin with, we point out that the usual gauge invariant 
perturbation theory cannot describe the fluctuations that we 
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actually observe. This is because we can observe only the 
fluctuations within the region causally connected to us. To 
discuss the so-called observable quantities in the framework 
of the gauge invariant perturbation, in general, it is necessary 
to fix the gauge in all region of the universe. However, in re- 
ality it is impossible for us to make observations imposing the 
gauge conditions in the region causally disconnected from us. 
Since we cannot specify the gauge conditions in the causally 
disconnected region, the gauge invariant variables that we 
usually consider as observables are undetermined. We need 
to be careful also in defining what are the observable fluctua- 
tions. We usually define the fluctuation by the deviation from 
the background value which is the spatial average over the 
whole universe. However, since we can observe only a finite 
volume of the universe, the fluctuations evaluated in such a 
way are inevitably influenced by the information contained in 
the unobservable region. In particular, in the chaotic inflation 
the longer wavelength mode has the larger amplitude of fluc- 
tuation lll|2|,|3l|4j], and therefore the value averaged over the 
whole universe is not even well-defined. In general, the devia- 
tion from the global average is much larger than the deviation 
from the local average, which leads to the over-estimation of 
the fluctuations due to the contribution from long wavelength 
fluctuations. 

In this paper, we show that, taking an appropriate gauge, we 
can compute the evolution of fluctuations which better corre- 
spond to what we actually observe. It is often the case to adapt 
the flat gauge or the comoving gauge in computing nonlinear 
quantum effects. Those are thought to be a way of complete 
gauge fixing. However, in § [II] we will explain that, even if we 
impose such gauge conditions in the observable finite region, 
the gauge conditions are not completely fixed. To remove the 
residual gauge degrees of freedom, we impose further gauge 
conditions. In doing so, we require also the gauge fixing con- 
ditions not to be affected by the influence from the causally 
disconnected region. The violation of causality due to care- 
less choice of variables, even if it is superficial such as pure 
gauge contributions, can lead to divergences in computation. 
In § [nfl we prove that IR corrections no longer diverge in the 
single field model, once we adopt an appropriate choice of 
variables with appropriate gauge conditions. We also show 
that the amplitude of perturbation does not grow secularly 
even if we send the initial time to the distant past unless very 
higher order perturbations are considered. In § |IV] we sum- 
marize our statement. 



II. A PRESCRIPTION TO SOLVE IR PROBLEM 



where M p i is the Planck mass. We perform the following 
change of variables 

= $/M pl , V{4>) = UW/Nfa, (2.2) 

to factorize M 2 ^ from the action as 

yp- r 

S = -f- J V=g [R - g^M» - 2V(<j))}d 4 x . (2.3) 

Hereafter we work with this rescaled non-dimensional field (j>. 
For simplicity, we assume that V((f>) and all of its higher order 
derivatives are at most 0(H 2 ), where H is the Hubble param- 
eter. This condition is satisfied in slow roll inflation 1 . Since 
the Planck mass is completely factored out, the equations of 
motion do not depend on it. The Planck mass appears only 
in the amplitude of quantum fluctuation. Namely, the typi- 
cal amplitude of fluctuation of $ is H, and hence that of <fi is 
H/M pl . 

In order to discuss the nonlinearity, it is convenient to use 
the ADM formalism, where the line element is expressed in 
terms of the lapse function N, the shift vector iV a , and the 
purely spatial metric h ab : 

ds 2 = -N 2 dt 2 + h ab (dx a + N a dt)(dx b + N b dt) . (2.4) 

Substituting this metric form, we can denote the action as 



Vh N (3) i? - 2NV{<j)) + — {E ab E ab - E 2 ) 



1 

N 



{4> - N a d a ^) 2 - Nh ab d a <t>d b 4> 



d 4 x 



where 



E ab = -{K b - D a N b - D b N a } 
E = h ab E a b . 



(2.5) 

(2.6) 
(2.7) 



In the ADM formalism, we can obtain the constraint equations 
easily by varying the action with respect to N and N a , which 
play the role of Lagrange multipliers. We obtain the Hamilto- 
nian constraint equation and the momentum constraint equa- 
tions as 



2V 

T-2i 



N~ (E E ab — E ) 
— N~*{cj) — N a d a (f)) 2 — h ab d a (, 

D a [N-\E a b - S a b E)} - N- l d b cj> 



%<t> = , (2.8) 
- N a d a <f>) = . 

(2.9) 



A. Setup of the problem 



Hereafter, neglecting the vector perturbation, we denote the 
shift vector as N a = d a x- In this paper we work in the flat 
gauge, defined by 



We first define the setup that we study in this paper. We 
consider the single field inflation model with the conventional 
kinetic term. The total action is given by 



2U{<i>)}d 4 x . (2.1) 



(2.10) 



1 This condition is not satisfied for small field inflation models. In that 
case we can relax the condition tod™ V((f>)/d<f> n = o(H 2 (M pl / H) n ~ 2 ) 
without changing the details of our arguments. 
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where a = e p is the background scale factor. Here we have 
also neglected the tensor perturbation, focusing only on the 
scalar perturbation, in which the IR divergence of our interest 
arises iH Si. 

In this gauge, using N, x an d the fluctuation of the scalar 
field ip, the total action is written as 



S 



1 



n=0 



+ N- 1 {-6p 2 + 4pA x 

+ (V Q V b X V a V & x-(Ax) 2 )} 
+ N-\4> + tp - V a X V Q ^) 2 - iV(V^) 21 



(2.11) 



and two constraint equations are 

n=0 

+ 4pA x + {V Q V b xV a V b x ~ (Ax) 2 } 
+ (0 + - V a xV a ^) 2 + iV 2 (V^) 2 = 0, 

(V a iV){2^ + (V a V bX ~ -5 Q b Ax)} 
-{V b <p)N (0 + 0- V a X V a ^) =0, 



(2.12) 



(2.13) 



where 



represents the three dimensional partial differentiation with 
respect to the proper length coordinates e p x and 

A = <5 o6 V V& . 

Spatial indices, a, are raised by S ab . This notation, 

which respects the proper distance, is convenient for the later 



discussions because it eliminates all the complicated scale fac- 
tor dependences from the action. 

The background quantities p and <fi satisfy 

■ip 2 = + v{4>) , 



4> + 3p <\> + = , 
where V$ = d^V ((f)). Expanding N, \ and ip as 



(2.14) 



N = 1 + <JiVi + -SN 2 + 
1 

X = Xi + H . 

1 



(2.15) 

we find that the first order constraint equations are written as 



V^tpx + 2VSN 1 + 2pAxi + <t> tp x = , (2.16) 

V a (2p 8N t - 4>(pi) = , (2.17) 
and the second order ones as 

+ 2VSN 2 + 2pAx 2 + 4> ^2 

+ 2pA X2 + V a V b xiV a V b xi - (Axi) 2 

- 20e-^V a xiV a (^ 1 + ipf + (V^) 2 = ,(2.18) 

V Q (2p 5N 2 -<j>ip 2 ) 

+ 2(V b <57V 1 )(V b V a xi-^ a Axi) 
-2V„^i(05iVi + 0i) =0. (2.19) 

Taking the variation of the action with respect to ip, we 
can derive the equation of motion for ip, which includes the 
Lagrange multipliers SN and x- F° r example, from the 
third order action, we can derive the equation of motion with 
quadratic interaction terms as follows, 



ip + 3ptp - Aip + V H <p - 4>A X + SNVj, - 3p<f>5N - d t {5N(j>) 

+~VM<f ~ V a (0 - j>5N)V a X - (0 - <j>6N)A x - pV a xV aV ? - d t {V a X V a v) 
-SpipSN ~d t (5Nip)~V a (5N\7 a ip) + V 4>( pipSN + 3p^5N 2 + d t (<p5N 2 ) = . (2.20) 

I 



Solving the constraint equations for the lapse function and 
shift vector at each order, we can express SN and x as func- 
tions of ip. Substituting these expressions into the original 
action ( 12. lit , we obtain the reduced action in the flat gauge 
written solely in terms of the dynamical degrees of freedom, 
(p. 



B. Tree-shaped graphs 

In this subsection, as a preparation for computing rt-point 
functions of <p(x), we consider an expansion of the Heisen- 
berg field tp(x) in terms of the interaction picture field <pi(x). 
When we compute n-point functions for a given initial state, 
we often use the closed time path formalism ||42[ l43i l44ll . in 



which n-point functions are perturbatively expanded by us- 
ing the four different types of Green functions: the Wightman 
function G + (x,x'), the Feynman propagator Gjr(x,x') and 
their complex conjugations G~(x,x') and Gd(x,x'). Here 
we shall adopt a different approach in which we take the full 
advantage of using the retarded (or advanced) Green function. 
In contrast to the above four Green functions, the retarded 
Green function 2 

G R (x,x') = iO{t - t')M^{G + (x,x') - G-(x,x')} , 

(2.21) 

is non-vanishing only when x is in the causal future of x'. 
In fact, when these two points are mutually space-like, the 
two field operators ip(x) and f{x') commute with each other, 
which leads to G + (x, x') = G~(x,x'). Since the retarded 
Green function Gr(x,x') has a finite non-vanishing support 
for fixed t and t', its three dimensional Fourier transform be- 
comes regular in the IR limit, while the other Green functions 
behave like fc~ 3 . (This IR behaviour leads to the scale invari- 
ant power spectrum, P(k) oc 1/fc 3 . We will discuss these 
issues in more detail later.) Hence, in order to prove the IR 
regularity in loop corrections to rt-point functions, it is conve- 
nient to use Gr(x, x') as much as possible. 

Let us denote the equation of motion for tp schematically as 

Lip. = -T[<p] , (2.22) 

where L is a second order differential operator corresponding 
to the linearized equation for ip (Eq. ( 12.371 )) and T stands for 
all the nonlinear interaction terms. We stress that this equation 
of motion does not depend on M v \ as anticipated. Using the 
retarded Green function Gr(x 1 x') that satisfies 

HG R {x, x') = -a,- 3 S 4 (x - x) , (2.23) 

we can solve Eq. ( 12.221 ) formally as 

ip(x) =(fi(x) + J d 4 x' G R (x,x')a 3 {t')T[ip](x'). (2.24) 

Here the factor a 3 originates from the background value of 
yj— g. Substituting this expression for <p{x) iteratively into 
T[tp] on the r.h.s., we obtain the Heisenberg field (f(x) ex- 
panded in terms of <pi(x) to any order using the retarded 
Green function Gr(x,x'). A diagrammatic illustration as 
given in Fig.|2]will be useful. In Fig.|2]we showed the proce- 
dure of expanding the Heisenberg operator when only a sim- 
ple three-point interaction is present, i.e. r[<p] = ■ Here, 
we represent the Heisenberg field, the interaction picture field, 
and the retarded Green function by a thick line, a thin line and 



4 

< +^c^ 




FIG. 2: Diagrammatic expression for the Heisenberg field expanded 
in terms of the interaction picture fields when only the three point in- 
teraction vertex is present. Here the Heisenberg field, the interaction 
picture field, and the retarded Green function are represented by a 
thick line, a thin line and a thin line with the index "R", respectively. 



a thin line associated with the index "R", respectively. Now it 
will be easy to understand that the Heisenberg field can be ex- 
pressed by a summation of tree-shaped graphs of this kind not 
only for this specific case, but also for any polynomial inter- 
action. Let us summarize the structure of tree-shaped graphs. 
Looking at a tree-shaped graph from left to right, it starts with 
a retarded Green function except for the first trivial graph that 
does not contain any vertex. All the retarded Green functions 
Gr(x, x') are followed by two or more <pi(x') or Gr(x' , x") 
with some integro-differential operators. All the interaction 
picture fields ipi (x) are located at the right most ends of the 
graphs. 

When we compute the expectation value for rt-point func- 
tions of the Heisenberg field, the interaction picture fields <pi 
that appear at the right ends of tree-shaped graphs are con- 
tracted with each other to make pairs. Then, when we eval- 
uate the expectation value, the pairs of ipj are replaced with 
Wightman functions, G + (x,x') or G~(x,x')(= G + (x',x)). 
As these propagators are IR singular (oc 1/fc 3 ) in contrast to 
Gr(x,x'), they are the possible origin of IR divergences in 
momentum integrations. 



C. Gauge degree of freedom in flat gauge 



2 Here G + and G~ are dimensionless propagator defined by G + (x, x') = 
(tpi(x)ipj(x')) and G~ (x, x') = (ipj(x')tp[(x)). Reflecting the overall 
factor Af 2 j in the action, these propagators are suppressed like 1/Af 2 j. As 
we use tne retarded Green function Gr to solve the equation of motion 
perturbatively, it is more convenient to define Gr not to dependent on 
M pl . Hence, M 2 j is multiplied in Eq. )2.2U . 



We consider the time evolution for the period [U, tf), where 
t f represents the final time at which we evaluate the field fluc- 
tuations. In this paper we assume that the universe is still in- 
flating at t = tf. Reflecting the fact that our observable region 
is bounded, we evaluate only the fluctuations within a finite 
region t/ , and we denote the causal past of this region t/ 
by 0. To exclude the effect from the unobservable part of the 
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universe, the evolution of tp in should be determined with- 
out any knowledge about the region outside 0. If T[ip] were 
written in terms of local functions of tp, Eq. ( 12.221 would de- 
termine the Heisenberg field tp{x) for x <G 0t f solely written 
in terms of the interaction picture fields tpi(x') with x' 6 0. 
However, we also need to solve equations of elliptic-type such 
as Eqs. dQ6b . dPTl ), d2~T8l ), and d2~T9l ). The solutions of 
these constraint equations, which determine the lapse func- 
tion and the shift vector, depend on the boundary conditions 
when a finite volume is assumed. Irrespective of the distance 
to the boundary, the boundary conditions immediately affect 
the solution owing to the non-hyperbolic nature of the equa- 
tions. 

In the linear order, these extra degrees of freedom appears 
as an arbitrary time-dependent integration constant. Indeed, 
we can solve the first order momentum constraint equation 
d2~T7T i as 

SN^x) = / x (i) + ^i(z) , (2.25) 
2p 

where an arbitrary function fi(t) was introduced as an in- 
tegration constant. Substituting this into the lowest order 
Hamiltonian constraint ( I2.16K we can solve it for \i to ob- 
tain 

where r is the proper spatial distance from the center defined 
by 

r 2 = e 2p x a x a . 

Since the last term in Eq. (I2.26l l proportional to r 2 cannot be 
expanded in terms of the spatial harmonics e lkx ), we do 
not have this residual gauge degree of freedom in the stan- 
dard cosmological perturbation scheme. Here, we do not care 
about the region outside 0. Then, the solution of Eq. ( 12.16b 
restricted to the region is not uniquely determined. Al- 
though we could have added more arbitrary harmonic func- 
tions (homogeneous solutions of the Poisson equation) with 
time-dependent coefficients to the above solution for xi> we 
neglected them for simplicity. 

The degree of freedom /i (t) introduced above corresponds 
to scale transformation: 

x a — ► x a = eP a[t) x a . (2.27) 

Such a scale transformation is compatible with the pertur- 
bative expansion only when our interest is concentrated on 
a finite region of spacetime. Once we consider an infinite 
volume, this transformation does not remain to be a small 
change of coordinates irrespective of the amplitude of a(t). 
Simultaneously, we apply the time coordinate transformation 
t — > t = t — a(t). Under this transformation, in the linear 
order, the spatial metric components are transformed to 

h ab (x) = e- 2 ^ h ab (x) 



= e 2 P^5 ab . (2.28) 

Thus, we find that this scale transformation keeps the flat 
gauge conditions that we imposed on the spatial metric (12. lot 
unchanged, and therefore it is in fact a residual gauge degree 
of freedom. 

Under the same coordinate transformation with the identi- 
fication 

ft 

h=a-^a, (2.29) 
2p 

we can easily confirm that the first order lapse function and 
the shift vector transform as given in Eqs. ( 12.251 and ( 12.261 >. 
Here we have explained only for the first order lapse func- 
tion and the shift vector, the corresponding degree of freedom 
also exists in the higher order. In this paper we focus on the 
flat gauge, but a similar discussion applies for the comoving 
gauge, too. 



D. Iteration scheme and local gauge conditions 

As is mentioned in § |U our final goal is to define finite ob- 
servable quantities in place of the naively divergent quantum 
correlation functions. We should note that in general, we can- 
not discuss observables in the gauge invariant manner by fix- 
ing the gauge completely over the whole universe. In this 
subsection, we show that imposing the boundary conditions 
unaffected by the information in the outside region, we can 
shut off the influence from the unobservable region of the uni- 
verse. (We refer to such a gauge as a local gauge, in which 
the causality is maintained also for the evolution of quantum 
Heisenberg field operators.) Once we choose the local gauge, 
we need not to care about the evolution outside the observable 
region as well as the gauge conditions there. 

Keeping the flat gauge conditions, we impose an additional 
local gauge condition: 

Wt<p(t) = jp J d 3 x W t (x) (p(t, x) = Q , (2.30) 

by using the degree of freedom fi (t) introduced in the preced- 
ing subsection and its higher order extension, where Wt(x) 
is a window function, which is unity in the finite region 
f = n Ef with a rapidly vanishing halo in the surrounding 
region, where S t means a t =const. hypersurface correspond- 
ing to the time t. For definiteness, we introduce 0' t D 0^ 
and define 0' as the causal past of 0' t . We require Wt(x) 
to vanish in the region outside 0'. In addition, Wt(x) is sup- 
posed to be a sufficiently smooth function so that an artificial 
UV contribution is not induced by a sharp cutoff. Lt, an ap- 
proximate radius of the region Oj, is defined such that the 
normalization condition 

W t l = 1 , 

is satisfied. 
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Roughly speaking, L t follows the radial null geodesic equa- 
tion. Hence, we have 



dt 

W) 



1 



Ltf + a{t)H{t) 



(2.31) 



For t > t c , we have L t ~ L tf , where t c is defined by 
a(t c )H{t c )L tf — 1. While, for t < t c , L t agrees with the 
comoving horizon radius at that time. (See Fig . |4] ) 

By construction, dp represents the deviation from the local 
average value in t . We associated "~" with the variables in 
this particular gauge, in order to clearly distinguish them from 
the variables for which the above additional gauge condition 
is not imposed. The difference between the variables with and 
without "~" is only in the boundary conditions. Hence, they 
obey the same differential equations, (12.171 i-( 12.201 i. 

Now we give a prescription to fix the arbitrary function 
fi(t) in Eqs. ( |2.25t and ( |2.26t as well as its higher order 
counterpart f n (t) (n — 2, 3, 4, • • • ) to satisfy the gauge con- 
dition (12.30b . For this purpose, we need to obtain a formal so- 
lution for dp. First, we consider the equations to fix the lapse 
functions. The higher order lapse functions are determined by 
the momentum constraint given in the form 



V a 5N n - —ip v 
\ 2p 



(n = l,2,3, 



(2.32) 



where the r.h.s. is a three vector at ?i-th order nonlinear terms 
expressed in terms of the lower order lapse functions, shift 
vectors, and (p. These equations do not have a solution in gen- 
eral since there are three equations with one variable. This 
situation happens because we have neglected the vector per- 
turbation. Hence, we consider only the scalar part of these 
equations, i.e. its divergence. This prescription is consistent 
with our neglecting the vector perturbation. The scalar part of 
Eq.(2.32) is formally solved as 



with 



SN n = SN n + /„, 



2p 



(2.33) 



The operation A -1 in Eq. ( I2.261 l is also to be defined so as 
to be completely determined by the local information in the 
neighborhood of t - We therefore define A -1 by 

where we have used the proper length coordinates X = efx 
. Similarly, the higher order shift vectors satisfy the Hamilto- 
nian constraint in the form 




where C n on the r.h.s. is a function expressed in terms of the 
lower order lapse functions, shift vector and (p. A formal so- 
lution for Xn is given by 



Xr, 



Xr, 



r 2 V 



fl I 



(2.35) 



with 



Xn 



-A~ 




->ri 



— A" 1 (v a E in) 
b V 



Next, we consider the equation of motion for dp, which is 
Eq. d2.20t with all perturbation variables replaced to the ones 
with Substituting the expressions for the lapse function 
( I2.331 l and the shift vector (12.35b into the equation of motion 
for dp truncated at the n-th order, we obtain an equation 



^Wn - 4>.fn + 



-y + 21/ J /„ 



-r„, 



where 



with 



L = dt + Zpdt - A 



A(t) 



~ 3p A 



(2.36) 



(2.37) 



and r„ on the r.h.s. of Eq. ( 12.361 1 represents all the remaining 
nonlinear terms expressed in terms of lower order terms in 
f(t) and dp. 

The equation for dp n is obtained by eliminating /„ from 
Eq. ( 12.361 ) by operating W t = 1 — W t as 



W t Hdp r , 



-W t T n [<p] 



(2.38) 



This equation alone is not sufficient to determine dp n because 
its homogeneous part is projected out. The homogeneous part 
of dp n is determined by the gauge condition W t <p n = 0. Prac- 
tically, dp n {x) is obtained by 



(2.39) 



(2.40) 



Lpn{x) = W t <p n (x) , 

where (p n satisfies 

&tp n (%) = -W t {x)T n [(p] . 

Here, for later convenience, we have inserted a window func- 
tion Wt(x) on the r.h.s. of Eq. d2.40| >, although it is possible 
to get the same conclusion without introducing this factor. As 
an effect of this inserted factor, thus obtained dp n (x) satisfies 
the field equation (12.361 1 only within the region 0. 

We found a way to obtain dp n before we know /„. Now we 
discuss how to fix /„. Operating Wt = J d 3 xWt(x) on 
Eq. ( 12.371 ), we obtain 




w t (r n + n<p n ), 



(2.41) 
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Using A(t), which satisfies the corresponding homogeneous 
equation 



4 = -^-2^ = 3^ 
A p (j) p 

we can solve Eq. d2.4U for /„ as 



24- 



(2.42) 



fn(t) 



A(t) 



dt' -±-Lw v (T n +SL(p n ). (2.43) 



Here we note that the r.h.s. is completely written in terms of 
the lower order perturbation variables and (p n , both of which 
are already given. 

From the above discussions we find that the lapse function, 
the shift vector and (p can be solved iteratively. Therefore all 
the higher order terms can be written in terms of tfii(x) = 

Wt<fi{x) with x £ 0' t .. In this sense, our prescription to 
find a solution of Heisenberg equations within guarantees 
approximate causality, avoiding influence from the outside of 
0'. We summarize our iteration scheme in Fig. [3] 

To summarize, we defined observable perturbations, which 
are not affected by the information in the region outside 0', 
by imposing an additional local gauge condition. Imposing 
appropriate boundary conditions in solving elliptic-type equa- 
tions that determine the lapse function and shift vector, we 
have shown that the local gauge condition that we require can 
be consistently imposed and the influence from the causally 
disconnected region is completely shut off in this gauge, in 
contrast to the traditional flat gauge. 



E. Quantization 

Even if we consider the perturbations in the local flat gauge, 
to quantize the fluctuation and to specify the initial state, we 
have to start with the ordinary flat gauge. This is because we 
need to give a complete set of mode functions on a Cauchy 
surface to specify the initial vacuum state and to constitute the 
Fock space. After specifying the initial state, we transform 
the perturbation variables into the local flat gauge, in which 



1 



( 



M2.39).(2.40)J^ 



f Eqs.(2.33), "\ 



> fn-f%n-f^n-\ 



the causality is maintained 3 . 

Using a set of mode functions {4>k(x) = Uk(t)e lkx }, which 
satisfy the linear perturbation equation 



= e- ikx L,j> k 

df + ipd t + e~ 2p k 2 



u k (t) , 
(2.44) 

we expand the globally defined interaction picture field as 

d 3 k (u k (t) ih< 



tpi(x) 



(2tt) 3 / 2 \ M p i 



afe + h.c. 



(2.45) 



Here the creation and annihilation operators, at and a k , sat- 
isfy the commutation relation 



a h ,a' k , =S 3 (k-k') 
The mode functions are normalized by 

Uk(t)u* k (t) - u k {t)u* k {t) = 



a3(t) 



(2.46) 



The initial vacuum state |0) is annihilated by the operation of 
any annihilation operator: 



a h \0) = 



for v fc 



We assume that the initial vacuum state is not so different from 
the adiabatic vacuum state at the initial time, especially for the 
long wavelength modes. On the initial surface, the Heisenberg 
operator corresponding to the scalar field fluctuation in the 
local flat gauge is related to that in the ordinary flat gauge as 
(p\{x) = W t ifi{x). 

III. IR REGULARITY 

In this section, we show that n-point functions for (p are IR 
regular for the most of inflation models. In § IIII Al we study 
the behavior of the mode function, especially focusing on the 
long wavelength limit. We show that the Wightman function 
is singular in the long wavelength limit like 1/fc 3 , while the re- 
tarded Green function is completely regular. In § IHIBI we will 
show that n-point functions calculated following our prescrip- 
tion are free from IR divergences in momentum integration 
when we do not care about secular growth of the amplitude of 
perturbation, which will be discussed in § IIII CI We will show 
that the secular growth does not occur unless very higher order 
perturbation is concerned. Secular growth means the increase 
of the amplitude of fluctuation in proportion to some power 



FIG. 3: Summary of the iteration scheme to obtain higher order per- 
turbation tp n . 



3 This gauge transformation on the initial time slice can be performed unam- 
biguously, because at the initial time we can safely neglect the non-linear 
interactions. Thus, we need not care about the ambiguity originating from 
the operator ordering. 
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of the e-folding number N in the slow-roll limit. Here, as we 
discuss more general setup in which H(t) monotonically de- 
creases, t is bounded by the condition H(t) < Hi <C M p \. 
Namely, the initial time is never sent to the infinite past where 
quantum gravity effects cannot be neglected. Therefore the 
e-holding number N in our setup is not infinitely long. In this 
sense, there arise no divergences from the time integration. 
Instead, our main concern is the dependence of the final result 
on the initial time ts. 



A. IR limit of mode functions and retarded Green function 

In this subsection we first discuss generic behavior of mode 
functions in the long wavelength limit. Using that general no- 
tion, we discuss the asymptotic behavior of the retarded Green 
function Gr(x,x') in the IR limit. 

To obtain the mode functions, Vk = ~*T"u k and the con- 

<p 

formal time coordinate rj = J dt/a are used. In terms of v k , 
Eq. fl2.44l ) becomes 



..(<*) 



+ k 2 v k ( V ) = , 



(3.1) 



where ' denotes a differentiation with respect to rj and 
z 2 = a 2 ((j)/p) 2 = —2a 2 p/p 2 . The normalization condition 
of mode functions ( 12.46b becomes 



z 2 (rj) 



(3.2) 



In the long wavelength limit we obtain two independent 
growing and decaying solutions as 



,(s) 



1 + fe 



(d) 

V, = — 
k 2 



" drf 

z 2 (r\') 
I r _drf_ 

Z 2 (l]') 



drj" z 2 ir\') + 



(3.3) 



Combining these two solutions, we can construct a mode 
function that satisfies the normalization condition ( 13. 2\ as 



1 



Vk 



W) vi9)+lC * ik) 



,(<*) 



with an arbitrary parameter c(k). 

To proceed further, let us consider a simple case in which 
the scale factor evolves as H = Hoa~ e , where e = —H/H 2 
is one of the standard slow roll parameters, and we assume 
that e is constant. Since the Hubble parameter should decay 
as a increases, e > is understood. As we are interested 
in the universe in an accelerated expansion phase, a oc a 1 ~ t 
should grow as a increases. Hence, e < 1 is also required. 
In this case, the original mode function u k is related to v k as 
Uk = — \/2ev k . The above two long wavelength solutions 
d3.31 > are reduced to 



,(fl) 



1 



2(1 



2 H 2 



2e(l 



)a 3 H 



(3.5) 



At the horizon crossing, where k w aH, the growing and 
decaying solutions should contribute to the positive frequency 
function v k to the same order if the initial quantum state is not 
very different from the adiabatic vacuum. Assuming that e is 
not very close to 1, this requirement determines the order of 
magnitude of c(k) as 



c(k) 



(3.6) 



First of all, from the above estimate of c(fc), we find that the 
leading order term in Uk in the long wavelength limit behaves 
like « H/V ek 3 . Hence, the Wightman function G + {x, x') = 
(ifi(x)(pi(x')) and its complex conjugation G~(x,x') have 
IR divergence. In fact, the Fourier transform of the Wightman 
function is given by 



(<pl(tM(t')) = u k (t)uUt'), 



(3.7) 



and it is 0(H 2 / (ek 3 )) in the long wavelength limit. 

The amplitude of oscillations of Uk changes approximately 
in proportion to 1/z oc l/a on sub-horizon scales, where 
k ^> aH. Hence, the amplitude of the positive and nega- 
tive frequency functions is enhanced for shorter wavelength 
modes compared with that in the long wavelength limit. How- 
ever, if such an enhancement causes problematic divergences, 
such divergences should be attributed to the issue of UV reg- 
ularization, which is not our main concern in this paper. On 
the other hand, the long wavelength limit of the decaying so- 
lution v/f \ grows faster than l/a as we decrease a. Hence, 

the absolute magnitude of the expression for in the long 
wavelength limit ( 13.51 ) gives an approximate upper bound on 

the true value of \vt \. 

In § IIHCI we will also use the expression for the retarded 
Green function Gr(x, x'). Formally, in terms of mode func- 
tions, we can give an expression for the retarded Green func- 
tion as 



G R (x, x') = -i6(t - t') 



d 3 k 
J^nj 3 



e lk < x - x '^R k (t,t') , 



(3.4) where 



(3.8) 



(3.9) 



Rki^t 1 ) = u k {t)u* k {t') - u* k {t)u k {t') . 
Substituting the expression ( 13.41 ). we obtain 

Rk(t( V ),t'( v ')) 

= -2ie (4 S) (rfV) " v k d Hv)vl g \v')) -(3.10) 

Then, using the expressions in Eq. ( 13. 51 ), we find that R k is 
regular in k without any singular behavior in the limit k — > 0. 

B. Momentum integration 

Now we are ready to discuss the IR regularity of n-point 
functions of (p(x). Our discussion is restricted to the case 
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excluding the slow roll limit e — > 0. In this limit, the 
background scalar field stays constant. Therefore we cannot 
choose Wtip = by a simple change of time coordinate. As 
a result, a singular behavior appears in Eq. ( 12.43b . In the fol- 
lowing discussion we do not care about the factor e in the final 
estimate of the order of magnitude, assuming that e is not ex- 
tremely small 4 D 

In this subsection we do not consider the secular growth of 
the amplitude of perturbation due to the integration for a long 
period of time. Namely, we consider the case that t, is not 
very distant past from t / . Therefore we do not care about the 
time integration. We defer this issue to the succeeding sub- 
section. Here we just consider the IR divergences originating 
from the momentum integration. We show that, if we follow 
the prescription described in § UH the amplitude of perturba- 
tion is IR regular without introducing any IR cutoff scale by 
hand. 

As (p{x) is composed of (p n (n — 1,2,3,- ••), we use 
the mathematical induction to show the regularity of all (p n . 
(p n (x) is, by definition, n-th order in the interaction picture 
field ipj. Formally, we define C[ip n ](x;p 1 ,- ■ ■ ,p n ) by ex- 
panding tp n (x) as 



n 



\J (2^)3/2"^ 



C[<p n ](x;Pi,-- ■ ,P n ) 

+ •••, (3.11) 



where we have suppressed terms containing creation opera- 
tors. The above expression is the result that we obtain after 
conducting all the integrations over the intermediate vertexes. 
The momenta {Pj} in the argument of C[0 ra ] are those associ- 
ated with the right most ends of the corresponding tree-shaped 
graph. 

What we will show below is the following properties of 

C[<Pn](x]Pi,-- ■ ,P n )- 

• It is a smooth function with respect to x in 0' for y pj = 
\pj | < a(t)A, where A is a momentum cutoff scale. 

• It vanishes in the long wavelength limit pj — ► 0. 

If C[ip n ] satisfies the properties mentioned above, one can 
easily show that n-point functions ((p(tf,Xx) ■ ■ ■ <p(tf,x n )) 
are free from IR divergences. When we take the expectation 
value of the product of (f>j (j < n) in the form of Eq. (13.1 lb . 
we consider all the possible ways of pairing with at,. 
Then, each pair of and at, is replaced with <5 3 (fc — fc'). One 
of the momentum integrations over k and k' is performed to 
obtain an expression in the form 



d 3 k 
(27rfc) : 



■C[<p ni ](xV, ■ ■ ■ ,fe,- • -)C[(p n2 ]{X2\- ■ ■ ,*!,••■)• 



4 It is well-known that T is suppressed by the slow-roll parameters. Hence, 
even in the limit e — » the integral in Eq. fl2.43t does not diverge, as long 
as the other slow-roll parameters scale in proportion to e. 



The remaining momentum integration does not have IR 
divergences owing to the second property of C[(p n ], i.e. 
limfe-,0 C[{p n ]{x; ■ ■ ■ , k, ■ ■ ■ ) = 0. 

For brevity, we denote a function which satisfies the above- 
mentioned two properties by an IR vanishing smooth function 
(IRVSF). In the following process of mathematical induction 
to show these properties, there is no operation on the momen- 
tum arguments in C[(p n ] . Not to confuse the readers, we stress 
that only the first argument, x, is relevant in the following dis- 
cussion. Our discussion in the rest of this subsection will pro- 
ceed mostly in the real space representation without switching 
to the Fourier space representation, because the finiteness of 
the volume 0' t is the clearer in the former representation. 

It will be obvious that IRVSFs satisfy the following proper- 
ties: 

Lemma If C\(x; {pj}) and C%(x',{qj}) are IRVSFs and 
there is no overlap between the list of momenta 
{p^ and {q 3 }, then V Q Ci(a;; {Pj}), x C x (x; {Pj}), 

Ci(x;{ Pj }), A^dfaipj}), Wtdfaipj}), 
J dtCifaipj}), and d(x; {p^) x d{x;{qj}) are 
all IRVSFs. 

To start the mathematical induction, one can easily check 

the first step that <pi(x) is an IRVSF. (p%(x) = Wt<pi{x) is 
expressed as 



where 



d 3 p 



(2tt) 3 / 2 



ip x _ Wt,—p 



U p (t) 



M, 



pi 



tip + {h.c.}, 



Wt.- P = J d 3 xe ip 



Wt(x). 



(3.12) 



(3.13) 



Here we note that Wt,o — J d 3 xWt(x) = L\. Hence, we 
have 



C[(fi](x,p) = 



W t ,_ p 



w. 



I A) 



P 3/2 u p (t) 



pi 



(3.14) 



This expression for C[(pi](x,p) is manifestly regular for the 
argument x. Inthelimitp — > 0, the factor [e lp ' x — W t ,- P /W ] 
vanishes. While the combination p 3 / 2 u p (t) is regular from 
the discussion in § IIII Al (See Eqs. ( 13.4b . ( 13.5b and ( 13.61 ).) 
Therefore C[<pi](x,p) vanishes in the limit p — > 0. Thus we 
find that C[<^i](x,p) is an IRVSF of 0(H/M p i) on super- 
horizon scales, p < a(t)H(t). 

The n-th order perturbation is obtained by 

£n = W t G R {W t >Y n ) 

= W t f dt 1 /dVa 3 (i')Gfl(x, a; ')^(x , )r„(a:'). 



(3.15) 

Wt' (x')T n (x') is constructed from lower order perturbations 
SNj, Xj, fj and (pj with j < n using the operations listed in 
the above Lemma. Furthermore, from Eqs. ( 12.331 ), ( 12.351 ) and 
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( I2.43l l. we find that SNj, Xj and fj are all constructed from 
(f>j by the operations listed there, too. Hence, C[Wi>r n ], the 
expansion coefficient of Wt'(x')T n (x') analogous to C[y3„] 
in Eq. ( 13.1 lb . is also an IRVSF. Since the expression of the 
retarded Green function ( 13.8b with Eq. ( 13.10b is regular in 
the IR limit, its Fourier transform Gr(x, x') should be reg- 
ular, too. (Regularity in UV is assumed to be guaranteed 
by an appropriate UV renormalization.) Since the integra- 
tion volume of x' is finite, the integral of a product of reg- 
ular functions Jd 3 x'a 3 (t')G R (x,x')Wt>(x')T n (x') should be 

finite, and hence it is IRVSF. Since the operation Wt preserves 
the properties of IRVSF, (p n = W t Gn(W t 'T n ) is also found 
to be IRVSF. 



C. Time integration 

In the preceding subsection we have shown that the ampli- 
tude of perturbation is regular as long as we do not care about 
the possibility of its secular growth. However, if we try to 
send the initial hypersurface to a very distant past, an- 
other significant amplification of the amplitude may arise. In 
this subsection we discuss this remaining issue, i.e. the initial 
time dependence of the amplitude of C[<p n ]. We will show 
that there is no significant secular growth in (p n for 

1 

n < n c — 1 , 

e 



and its amplitude is bounded by 



c[q> n \ < o{A n ), 



(3.16) 



where 



H 



(ciiHiLt) 



M„ 



for n < n c , 
for n > n c , 



Time integration appears not only in Eq. ( 13.21b but also 
in Eq. ( 12.43b . Both contain the interaction vertex r„. The 
interaction vertexes and the retarded Green function Gr do 
not contain M v \ because the the factor Mpj is completely fac- 
tored out in the action. Hence, all the dimensional coefficients 
whose mass dimension is one are 0(H). Owing to the as- 
sumption of induction, we have ipj — 0(Aj) for j < n. 
Thus, based on dimensional analysis, the order of magnitude 
of C[r„] is estimated as 



C[T„] = H 2 Q (£ r . 



(3.17) 



with 



H 



max- 



H 



A/„ 



,(a t H,L t )^- 



for n < n c , 
n-l 



Mr, 



(oiHiLt)* 



Hi 



M„ 



') 

^ , for n > n c , 



Here we have used 



C[f j ]=0('B j ) , 



(3.18) 



for j < n, which will be proven immediately below. 

To derive this rough estimate of the order of magnitude, we 
use the simple model introduced in § IIII Al again. Assuming 
that H oc a~ e with < e < 1 as before, we read the time 
integration in Eq. ( 12.43b as 



fj 



1 



< 



da(f) a 3 (t') 
2ea 3 H J a(t') H{?) 



(3.19) 



where in the last inequality we have assumed that the integral 
is dominated by the later epoch, t' w t. For j < n c , this is 
always the case. Then, it will be obvious that the condition 
(|3.18|) is satisfied. For j > n c , a similar argument holds when 
the integral is dominated by the later epoch. However, there is 
also a possibility that the integral is dominated by the earlier 
epoch. In this case we have 



fj 



< 



a(tfH(t) 



O 



H, 



1 1 



A I, 



pi 



(3.20) 



(Notice that aiH t L u ps 1. ) Since a 3 Hi/a 3 (t)H(t) < 
(ciiHiLt) 3 < (aiHiLt) 2 , the condition (|3. 18|) is satisfied in 
this case, too. 

Now we turn to the time integration in Eq. ( 13.15b . which 
can be expressed, using the Fourier component (Wt'T n )k = 



/ d 3 x' e- ik x 'W t> (x')T n (x'), 



G R (Wt>T n ) = 



as 



dt' 



d 3 k 



xa 3 (t')R k (t,t')(W t ,r n ) k . (3.21) 



(We have introduced Wt' in Eq. ( 12.40b in order to make 
the Fourier component (Wt'^njk well-defined here.) Since 
(Wt'T n )k is a regular function whose non-vanishing support 
is limited to a finite region, its Fourier coefficient is also reg- 
ular as a function of k. When we consider a fixed value of k, 
the time integration should be truncated at tk defined by 

k = a(t k )H(t k ), 

due to the U V cutoff 5 . Thus the relevant modes for the inte- 
gration over a long period of time are concentrated on small 



5 As we have mentioned in § 1111 Al there is an enhancement of the amplitude 
of tpi(x) = Wt<Pi{x) for sub-horizon modes. However, the momentum 
integration including T n should be dominated by the modes near the hori- 
zon scale or the modes with a longer wavelength. If the contributions from 
the shorter wavelength modes dominated, the results of computation would 
depend on the UV cutoff scale A. Then, some factors of H in the above 
estimate of the order of magnitude would be replaced with A. However, the 
appearance of A in the final results means that the UV renormalization has 
not been properly done. If the UV renormalization is appropriately con- 
ducted, the counter terms should cancel the contributions which increase 
toward the shorter wavelength modes so that the cutoff scale A do not ap- 
pear in final results. Then, the contributions from the sub-horizon scales 
do not affect the order of magnitude of <p n . This means that, owing to an 
appropriate UV renormalization, we can safely assume that the effective 
UV cutoff momentum scale is as small as H. 
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k limit. (In this sense, the problem of initial time dependence 
(or secular growth) is a kind of IR divergence problem.) Since 
the inequality k < a(tk)H(tk) < a(t)H(t) holds for the rel- 
evant modes in the inflating universe, we can assume that all 
the modes in the momentum integration are on super-horizon 
scales at t. Then, one can use the long wavelength expansion 
for Vk(t) in Eq. (13.10b . For v^'(t') the expression in the long 
wavelength expansion is not a good approximation. However, 
as we have seen in § IIII Al the leading order expression for 



,<<') 



(t 1 ) in the long wavelength limit can be used as an esti- 



mate of the upper bound of its magnitude. Thus we find 

1 



|ifc(t,Ol«|2e«£V)l 



< 



a 3 (t')H(t>) 



(3.22) 



Using this expression for the retarded Green function, 
Eq. (13. 2U is estimated as 



\{G R {W t ,T n )) k \ < 



a{t) Mf) 

a(t k ) 0(f) 



(w t ,r n ) k 



H 2 (t') 



(3.23) 



Using the fact that the amplitude of the Fourier coefficient 
{Wt'T n )k is bounded by the amplitude of T„(x') multiplied 
by the volume of the window function L\, . 



C[G R (W t 'T n )] k 



< 



a(t) 



da(t') , 



a(t k ) 0(f ) 



To proceed further, we divide the area of the above integra- 
tion in two dimensional space of (k, a(t')) into two regions; 
(i) k > L^ 1 and (ii) k < L^ 1 , as shown in Fig. [4] We dis- 
criminate the region (ii) in which the operation of Wt results 
in an additional suppression of amplitude from the region (i) 
in which it does not. 

Let us consider first the region (i). In this case, as is obvious 
from Fig. [4] the time integration is restricted to t' > t c . Hence, 
we have Lf ~ L t ~ L tf . Furthermore, since t c is not so 
far from tf, we can approximate H(t') by H(t). Hence, for 

logo 



\f\ct n . 


< 




log a(t c ) 


(ii) 


/ 

/ 

s 




/[ ~ > ] 






/ -iogLt f 



FIG. 4: The dark grey region represents the region (i). The light 
grey region represents the region (ii). These two regions are divided 
by the solid curve which shows the scale of the causally connected 
region, i.e., k — L^ 1 . The dashed line is the horizon scale, i.e., 
k = a(t)H(t), which corresponds to the effective UV cutoff scale. 



n < n c we obtain 



C[G R (W v T n )} h 



< 



a(t k ) 



a(t') 
H(t) 



L 



' V Mpi J 



pi 



log{a{t)/a(t c )) 



Since \og(a(t) / a(t c )) cannot be a large number, this inequal- 
ity means that C[G R (W t 'T n )} = 0{[H(t)/M p i] n ). A paral- 
lel argument holds for n > n c , too. 

Next, we consider the region (ii). To consider this region, 
it is essential to take into account the operation of Wt- Using 
the relation ( 13.24b . we have 



C[W t G R (W t ,T n )] 
da(t') 



< 



a(t>) 

ik-x 



d 3 k 
W t ,o . 



)4 


\H(t')] 




Mpi . 



(3.25) 



Here we have changed the order of integrations. 

As we consider the region k < L^ 1 , Wt — k/Wt,o will t> e 
expanded as Wt,-k/Wt,o = 1 + 0((L t k) 2 ). Therefore the 

factor W t e lk x — e lk x — W t ,-k/W t ,o is approximated by 
i k x . Then, performing the momentum integration, we obtain 



C[W t G R (Wt'T n ) 



dag) j. 
a(t>) W 



H(t>) 



(3.26) 



where we have used \x\ < L t . For small a(t'), L t , x w 
a(t')H(t'). Thus we find that the integrand L^\H/M p {) n is 
proportional to (aH)H n cx a 1_ (" +1 ) c . Forn < n c = e _1 — 1, 
the integration is dominated by the later epoch. Hence, we 

have an estimate C[W t G R {W V T n )] = O ([H(t)/M p i] n ). For 
n > n c = er 1 — 1, there are terms whose order of magnitude 
is bounded by 



or 



(atHiLt) 



(aiHiLt) 



Hi 



M, 



pi 



da(t') H(t') 



o(f) Mj 



pi 



Hi 



pi 



da(t') 



(aiHiLt'), 



besides the terms that are estimated as in Eq. ( 13.261 ). In 
all cases the time integration is dominated by the earlier 

epoch. Thus, we have an estimate C\W t G R (W t 'T n )\ = 
O ((aiHiLt) [Hi/M p i] n ) . We find that the initial time de- 
pendence remains in C[y> n ] for n > n c , but it has at least 
one suppression factor (aiHiLt) associated. To conclude, we 
have shown that the condition ( 13.16t is satisfied in all cases. 

Before closing this section, we would like to stress the im- 
portance of the factor Wte x in Eq.( |3.25l l, which is absent 
in the standard treatment. This factor is the origin of the factor 
L t /L t i in Eq. ( 13.261 ). If it were not for this factor, this integral 
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would be dominated by the earlier epoch for any n. (In the 
slow-roll limit e — ► 0, the integral would be proportional to 
the e-holding number N = log a(t)/a(ti).) Hence, the ini- 
tial time dependence appears in the n-point functions even for 
small n at the lowest tree-level order. 

Even if we follow our improved prescription, the contri- 
bution from (p n with n > n c carries the dependence on the 
artificial choice of the initial time ij. Physical origin of this 
dependence on t{ is clear. This dominance of the contribution 
from the earlier epoch originates simply from larger ampli- 
tude of fluctuation due to larger H. Even if the propagation of 
fluctuation from far past is suppressed, the source T n rapidly 
increases toward the past for large n. Therefore we suspect 
that this initial time dependence might be really physical, al- 
though it appears only when we consider sufficiently higher 
order perturbations. However, as we have not used all the 
residual gauge degrees of freedom, there might be a better 
prescription for the gauge fixing in which the critical order n c 
is larger. 



IV. CONCLUSION 

As the possibility of detecting nonlinearities in the primor- 
dial perturbations of the universe is increasing, it becomes 
more important to understand the issue of IR divergences in 
the computation of primordial perturbations and to predict 
their finite amplitude that we actually observe B45I1 . In this 
paper, we pointed out that the standard prescription of the cos- 
mological perturbation theory contains residual gauge degrees 
of freedom if the gauge conditions are imposed only locally 
within our observable universe, and that it is important to fix 
these gauge degrees of freedom to remove IR divergences. In 
order to fix the residual gauge degrees of freedom, taking the 
boundary conditions which shut off the influence from the un- 
observable region, we proposed the use of local gauge fixing 
conditions. 

When we have an equation of elliptic type, the boundary 
conditions are not arbitrary in general. If we change the 
boundary conditions for an elliptic type equation, we obtain 
a different solution. However, here the elliptic type equa- 
tions appear only for determining the lapse function and the 
shift vector. The boundary conditions in solving the ellip- 
tic type equations are not specified from the fiat gauge con- 
dition alone. A different choice of the boundary conditions 
corresponds to a different way of fixing the residual gauge de- 
grees of freedom. Our choice of local gauge conditions is not 
unique, but it completely fixes the gauge in without using 
any information outside 0'. 

It is true that the n-point functions calculated in the present 
manner depend on the choice of fixing the residual gauge. 
Making use of the transfer functions, any real observables like 
the angular power spectrum of the CMB sky map can be de- 
scribed in terms of these ri-point functions for the primordial 
perturbations H46I1 . For the single field inflation model, we 
have shown that the amplitude of our primordial perturbations 
is free from IR divergences (unless the Hubble parameter at 
the initial time is well below the Planck scale). Then, the real 



observables should be also IR regular. We also pointed out the 
possibility that the terms which depend on the initial time may 
dominate in higher order perturbations above a critical order. 

At the end of this paper, let us comment on the case in 
which more than one fields participate in IR divergences. In 
our proof of the absence of IR divergences we used the gauge 
in which the local average of the inflaton field does not fluctu- 
ate using one of the residual gauge degrees of freedom men- 
tioned above. This adjustment of the average value is possible 
only for one field. When plural fields have scale invariant or 
even redder spectra, therefore our prescription presented here 
is not enough to regularize IR divergences. This claim is on 
the same line with the ar gum ent given by G. Geshnizjani and 
R. Brandenberger in 1471 14811 . Discussing the backreaction on 
the background expansion rate due to classical fluctuations, 
they showed that the observable expansion rate does not suf- 
fer from cumulative backreaction in single-component mod- 
els, while it does in multi-component models. 

Thus, when plural fields are concerned with IR divergences, 
we need more careful discussion about what we actually ob- 
serve. When we consider the eternal inflation scenario, the 
wave function of the universe is infinitely spread in the field 
space, and the expectation values of field fluctuations will di- 
verge. We think that these divergences due to the fields other 
than inflaton are physical. However, in the actual observa- 
tion of the universe we will not see any divergences. The key 
idea will be that what we compute as the correlation func- 
tions in field theory are different from what we really observe. 
We think that in this case it is essential to take into account 
the decoherence effects in order to remove these IR diver- 
gences. Deferring the detailed explanation to the succeeding 
paper [49], we describe here our basic idea how to handle the 
divergences in the multi-field case briefly. We focus on the 
field whose IR corrections still diverge even after the local 
gauge fixing. We denote it by <pm. The adiabatic vacuum 
state can be decomposed into a superposition of wave pack- 
ets which have a peak at a certain value of the local average 
Wt<pin(Tf). As the universe evolves, the wave packets lose 
correlation to each other. Through this so-called decoherence 
process, the coherent superposition of the wave packets starts 
to behave as a statistical ensemble of many different worlds, 
where each world means the universe described by a deco- 
hered wave packet lf50l EH l52ll . Our observed world is just 
a representative one expressed by a wave packet randomly 
chosen from various possibilities. Once one wave packet is 
selected after the decoherence process, the evolution of our 
world will not be affected by the other parallel worlds. How- 
ever, the initial vacuum state does include the contributions 
from all the wave packets. This implies that a naive computa- 
tion of n-point functions is contaminated by the contribution 
from the other worlds uncorrelated to ours, which is the origin 
of the divergences. 

Recently the stochastic approach 01 III |M III III III III 
has been employed in order to solve the IR divergence prob- 
lem ll60ll6Tll62ll . This is in harmony with our claim. However, 
it is hard to deny the spiteful suspicion that the reason why the 
problem of IR divergence does not appear in the stochastic ap- 
proach might be simply because quantum fluctuations in the 
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IR limit are neglected by hand. Therefore, in our succeeding 
paper, we describe the decoherence effect without relying on 
the stochastic approach, and discuss the regularity of the IR 
corrections. 

In contrast, when we consider the case in which only 
a single field is responsible for IR divergences, using the 
residual gauge degrees of freedom, we can adjust the local 
average value of the field not to fluctuate. Then, as we have 
shown in this paper, we need not to pick-up one decohered 
wave packet from the superposition of infinitely many wave 
packets. 
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